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We extend the formalism of weak interaction processes, obtaining new expressions for the tran- 
sition rates, which greatly facilitate numerical calculations, both for neutrino-nucleus reactions and 
muon capture. Explicit violation of CVC hypothesis by the Coulomb field, as well as development 
of a sum rule approach for the inclusive cross sections have been worked out. We have done a 
thorough study of exclusive (ground state) properties of ^'^B and ^'^N within the projected quasi- 
particle random phase approximation (PQRPA). Good agreement with experimental data achieved 
in this way put in evidence the limitations of standard RPA and the QRPA models, which come 
from the inability of the RPA in opening the P3/2 shell, and from the non-conservation of the num- 
ber of particles in the QRPA. The inclusive neutrino/antineutrino {i^/i') reactions ^^C(!^, e^)^^N 
and "C(i>, e+)"B are calculated within both the PQRPA, and the relativistic QRPA (RQRPA). 
It is found that the magnitudes of the resulting cross-sections: i) are close to the sum-rule limit 
at low energy, but significantly smaller than this limit at high energies both for v and D, ii) they 
steadily increase when the size of the configuration space is augmented, and particulary for v/i^ 
energies > 200 MeV, and iii) converge for sufficiently large configuration space and final state spin. 
The quasi-elastic ^^C(!/, /^~)^^N cross section recently measured in the MiniBooNE experiment is 
briefly discussed. We study the decomposition of the inclusive cross-section based on the degree of 
forbiddenness of different multipoles. A few words are dedicated to the v jv-^'^Q, charge-exchange 
reactions related with astrophysical applications. 

PACS numbers: 23.40.-s, 25.30.Pt, 26.50.-|-x 
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I. INTRODUCTION 

The massiveness of neutrinos and the related oscil- 
lations are strongly sustained by many experimental 
works involving atmospheric, solar, reactor and accel- 
erator neutrinos [lH7[. The subsequent experimen- 
tal goal is to determine precisely the various parame- 
ters of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
neutrino mass matrix, absolute masses of different fla- 
vors of neutrinos, CP violation in neutrino sector, etc. 
To address these problems several analyses of neutrino 
oscillation data are presently going on. At the same 
time, several experiments presently collect data, and oth- 
ers are planned. Accelerator experiments, experiments 
with neutrinos from i/- factories, /3-beams, etc., are also 
planned and designed, as well as some experiments with 
natural i^-sources like solar neutrinos, atmospheric neu- 
trinos, or antineutrinos from nuclear reactors. 

The neutrino-nucleus scattering on ^^C is important 
because this nucleus is a component of many liquid scin- 
tillator detectors. Experiments such as LSND [l|, Qj 
KARMEN [i, Q, and LAMPF 0, [U have used ^^c 
to search for neutrino oscillations, and for measuring 
neutrino-nucleus cross sections. Present atmospheric and 
accelerator based neutrino oscillation experiments also 
involve ^^C, and operate at neutrino energies Ey ^ 1 
GeV in order to access the relevant regions of oscillation 
parameter space. This is the case of the SciBar detec- 



tor [I4I, where the molecule Cgi^s is involved, and the 
MiniBooNE detector [TjI, which uses the light mineral 
oil containing the molecule Cif2- The ^^C target will be 
used in several planned experiments, such as the spalla- 
tion neutron source (SNS) at Oak Ridge National Lab- 
oratory (ORNL) [ll, and the LVD (Large Volume De- 
tector) experiment [1^ , developed by the INFN in Gran 
Sasso. 

For the planned experimental searches of supernovae 
neutrino signals, which involve ^^C as scintillator liquid 
detector, the precise knowledge of neutrino cross sections 
of ^^N and "'^^B ground-states, i.e., of a^-(Ey,\\), and 
crp+ [Ej-j, 1 ^) is very important. In fact, in the LVD ex- 
periment [l5[ the number of events detected during the 
supernova explosion are estimated by convoluting the 
neutrino supernova flux with: i) the interaction cross 
sections, ii) the efficiency of the detector, and iii) the 
number of target nuclei. For the carbon content of the 
LVD detector have been used so far a^~(Ey,Y\^). and 
(Te+(i?i;, 1^), as obtained from the Elementary Particle 
Treatment (EPT) [l^. Moreover, as an update of the 
LVD experiment related to supernovae neutrinos detec- 
tion (where ""^^C will also be employed), there is ongoing 
design study concerning large size scintillator detectors, 
called LAGUNA, where a 50 kt scintillator LENA is be- 
ing considered [l3 • 

On the other hand, as the ^^C nucleus forms one of the 
onion-like shells of a large star before collapse, it is also 



important for astrophysics studies. Concomitantly, sev- 
eral authors [1^, [l8M23| have recently stressed the impor- 
tance of measuring supernova neutrino oscillations. They 
claim that a supernova explosion represents a unique sce- 
nario for further study of the PMNS matrix. The cor- 
responding neutrinos, which carry all flavors were ob- 
served in only one occasion (SN1987A), have an energy 
Ev < 100 MeV [11], and are also studied through the 
interactions with carbon nuclei in the liquid scintillator. 

Thus, the main interest in the neutrino/antineutrino- 
"'^^C charge- exchange cross sections comes from the neu- 
trino oscillations, and precise knowledge of the cross sec- 
tions in the neutrino energies going from a few McVs up 
to a few GeVs is required. Up to quite recently the only 
available experimental information on reactions was that 
for the flux-averaged cross-sections: i) ^^C(i'e,e~)^^N 
in the DAR region: E^ < 60 MeV [H-ijI, and ii) 

MeV [2i-i 

mental programs at MiniBooNE [311, K2K [32], and Sci- 
BooNE [3^ yield results on the {v^^^'^G) cross section 
for 0.4 GeV < E^^ < 1.7 GeV. It is weU known that 
for E,y larger than a few hundreds MeV's, besides the 
quasi-elastic (QE) channel, many inelastic channels are 
open and pion production becomes important. In fact, 
there have been quite active experimental efforts to inves- 
tigate neutrino-induced coherent single-pion production 
in the A-excitation region of ^^C. Starting approximately 
at the threshold coming from the pion and charged lep- 
ton masses (m^ and m^), the n + £ production cross sec- 
tion steadily increases with the neutrino energy becoming 
larger than the quasi-elastic one for Ei,^l.5 GeV [3ll - l33l | . 

From the theoretical side there have been great ef- 
forts to understand the nuclear structure within the 
triad {^^B,^^C,^^N}. In the seminal work of O'Connell, 
DoncUy, and Walecka [3J] a unified analysis of electro- 



i^N in the DIF region: 127 MeV < E^^ < 300 
In last few years, however, several experi- 



magnetic, and semileptonic weak interactions was pre- 
sented. To describe the nuclear dynamics they have used 
the particle-hole Tamm-Dancof Approximation (TDA) 
within a very small single-particle space ^ {S2 = {lsi/21 
lP3/2, lPi/2, lc^5/2, 2S1/2}) [11]. To achieve agreement 
with experiments for the /3^-decays, and /^-capture they 
were forced to use an overall reduction factor ^^ of the 
order of 4 (2) for even (odd) parity states. They have 
also pointed out that this factor would become totally 
unnecessary with use of a better nuclear model able to 
open the 1^3/2 shell. 

Rather thorough comparisons of 2sld and 2plf shell- 
model predictions with measured allowed /3-decay rates 
have yielded a simple, phenomcnological effective axial 
coupli ng q ^ = 1 that should be used rather than the bare 
value |36l439| . This observation is the basis for many 



^ From now on a singic-particle (s.p.) space that includes all or- 
bitals within A'^ harmonic oscillator (HO) shells will be labeled 
as space Spf. 



nuclear model estimates of the Gamow- Teller (GT) re- 
sponse that governs allowed neutrino cross sections. In 
Ref. [3J] g^ = 1.23 was used based on a study of neutron 
/3-decay, and, as the analyzed processes were dominantly 
of the axial- vector type, the use oi g^ = 1 would have di- 
minished the reduction factors ^^ in an appreciable way. 

In the Random Phase Approximation (RPA), besides 
the TDA forward-going amplitudes, the backward-going 
amplitudes are present as well. However, these additional 
RPA amplitudes did not help to open the IP3/2 shell 
in the continuum RPA (CRPA) calculations of Kolbe, 
Langanke, and Krewald [40]. Thus, as in the case of the 
TDA used in Ref. [3J] , to get agreement with data for the 
ground state triplet T = 1 (/3^-decays, /x-capture, and 
the exclusive ^'^C{ve,e~)^'^N reaction) their calculations 
were rescaled by a factor = 4. 

The main aim of the CRPA is to describe appropri- 
ately not only the bound states but also the virtual 
(quasi-bound), resonant, and continuum states, which 
are treated as bound states in the RPA. However, this 
superiority has not been evidenced so far in numerical 
calculations. For instance, in the case of ^-capture rates 
in ^^N the two methods agree with each other quite well 
for the 0^ and 1~ states, while the RPA result is pre- 
ferred for the 2^ state [4ll |. 

To open the IP3/2 shell one has to introduce pair- 
ing correlations. This is done within the Shell Model 
(SM) (42h44| . which reproduces quite well both i) 
the experimental flux-averaged exclusive, and inclusive 
cross sections for the ^^Cive^c-y^N DAR (H-H^, and 
i2C(i/p,/i-)^2j^ DIF in reactions, and h) the ^i~ + 
^^C — >■ f^ -I- ^^B muon-capture modes [45l - l47l |. 

The quasiparticlc RPA (QRPA) also opens the IP3/2 
shell by means of the pairing interaction. However, it 
fails as well in accounting for the exclusive processes to 
the isospin triplet T = 1 in ^^C, because a new problem 
emerges, as first observed by Volpe et al. [43|. They 
noted that within the QRPA the lowest state in ^^N 
irremediably turned out not to be the most collective 
one. Later it was shown (48l - [5(j | that: 1) the origin 
of this difficulty arises from the degeneracy among the 
four lowest proton-neutron two-quasiparticle (2qp) states 

|lPl/2lP3/2), |lP3/2lP3/2), \^Pl/2^Pl/2) and | lp3/2lPl/2) , 

which, in turn, comes from the fact that for N = Z = 6 
the quasiparticle energies Eip-^^^ and Eip^^^ are very close 
to each other, and 2) it is imperative to use the projected 
QRPA (PQRPA) for a physically sound description of 
the weak processes among the ground states of the triad 
{12b, 12c, i2n} [H-i^; see Figs. 2 and 3 in Ref. ^. 

In summary, neither the CRPA nor the QRPA are the 
appropriate nuclear models to describe the "fine struc- 
ture" of exclusive charge- exchange processes around ^^C, 
and they only can be used for global inclusive descrip- 
tions. Of course, the same is valid for the relativistic 
RPA (RQRPA) that has recently been applied with suc- 
cess in calculations of inclusive charged-currcnt neutrino- 
nucleus reactions in ^^C, ^^O, ^^Fe, and ^ospb [H, and 
total muon capture rates on a large set of nuclei from ^^C 



to Pu [52] ■ The continuum QRPA (CQRPA) would section and the angular distribution of the outgoing elec- 



have to be superior to the QRPA for the same reasons 
that the CRPA would have to be better than the RPA. 
Nevertheless, neither this superiority has been put in ev- 
idence by numerical calculations |53l . |5J|. Finally, it is 
clear that the nuclear structure descriptions inspired on 
the Relativistic Fermi Gas Model (RFGM) [HsI - fSTt . which 
do not involve multipole expansions, should only be used 
for inclusive quantities. 

When the effects due to resonant and continuum states 
are considered, as it is done within the CRPA and 
CQRPA, the spreading in strength of the hole states in 
the inner shells should also be taken into account for the 
sake of consistency. In fact, a single-particle state j that 
is deeply bound in the parent nucleus, after a weak inter- 
acting process can become a highly excited hole-state j^^ 
in the continuum of the residual nucleus. There it is sud- 
denly mixed with more complicated configurations (2hlp, 
3h2p, . . . excitations, collective states, and so on) spread- 
ing its strength in a relatively wide energy interval [58|^. 
This happens, for instance, with the lsi/2 orbital in ^^C, 
that is separated from the 1^3/2 state by approximately 
23 MeV, which is enough to break the 12 particle system, 
where the energy of the last excited state amounts to 11.5 
MeV in ^^n, and 16.5 McV in ^^B channels. Although 
the detailed structure and fragmentation of hole states 
are still not well known, the exclusive knockout reactions 
provide a wealth of information on the structure of single- 
nucleon states of nuclei. Excitation energies and widths 
of proton-hole states were systematically measured with 
quasifree (p, 2p) and (e, e'p) reactions, which revealed 
the existence of inner orbital shells in nuclei [59-1631 . 



In the TDA calculation of Ref. [3J] the 5*2 space has 
been used, which extends only from 13.77 MeV up to 
30.05 MeV, embracing, respectively, 1, 2, 2, 1, and 1 neg- 
ative parity states J'^ = 0~,1~,2~,3^, and 4^, and 
1, 2, 2, and 1 positive parity states J'^ = 0+, 1+, 2+, and 
3+. With such small configuration space, the neutrino 
cross sections ae{Ei,), and afj_{E^) have been evaluated 
up to a neutrino energy E^, of 0.6 GeV, and extrapolated 
up to 20 GeV. In recent years, however, large configu- 
ration spaces have been used in the evaluation of QE 
cross sections for E^ ^ 1 GeV. For instance, Amaro 
et al. [631 have employed the single-particle SM (TDA 
without the residual interaction) in a semirelativistic de- 
scription of quasielastic neutrino reactions (v^,^") on 
^^C going up to El, = 1.5 GeV, and including multi- 
poles J'^ < 47*. Good agreement with the RFGM was 
obtained for several choices of kinematics of interest for 
the ongoing neutrino oscillation experiments. Kolbe et 
al. [69| have also achieved an excellent agreement between 
the RFGM and the CRPA calculations of the total cross 



^ One should keep in mind that the mean life of ^^N and ^^B arc, 
respectively, 11.0 and 20.2 ms, while strong interaction times are 
of the order of 10~^^ s. 



trons in ^^0{i'e,e)X for E^, < 0.5 GeV. They have con- 
sidered states up to J^ ~ 9* only, and didn't specify the 
configuration space used. Moreover, Valverde et al. [53| 
have analyzed the theoretical uncertainties of the RFGM 
developed in [56[ for the (i/e,e~), and (t^^,/i~) cross sec- 
tions in i^C, i6n and ^OCa for E^ < 0.5 GeV. The work 
of Kim et al. [70| should also be mentioned where were 
studied the effects of strangeness on the (j/^,/i~) and 
{i'f^,li~^) cross sections in ^^C for incident energies be- 
tween 0.5 MeV and 1.0 GeV, within the framework of a 
relativistic single-particle model. Quite recently, Butke- 
vich [7l| has also studied the scattering of muon neu- 
trino on carbon targets for neutrino energies up to 2.8 
GeV within a relativistic shell-model approach without 
specifying the model space. 

For relatively large neutrino energies {E^^ ^ niTr , and 
-Bi/j, ^niTT+m^) to the above-mentioned QE cross sections 
should be added the pion production cross section, as 
done, for instance, in Refs. [72, uM- One should also 
note that ae{Ei,^) and (t^(_Ej/ J coincide with each other 
asymptotically. This is clearly put in evidence in the 
Extreme Relativistic Limit (ERL) where \\ii\/Ei — > 1, 
and the neutrino-nucleus cross sections depend on m^ 
only trough the threshold energy, as can be seen from 
the Appendix [C[ of the present work. The Figure 4 from 
Ref. [TJI is also illustrative in this respect. 

Therefore, we focus our attention only on the quasi- 
elastic cross section af,{E,^^) since, at muon-neutrino en- 
ergies involved in the MiniBooNE experiment [IjI, it is 
equal to ct^{Ei, ) for all practical purposes. 

One of main objectives in the present study is to ana- 
lyze the effect of the size of the configuration space up to 
neutrino energies of several hundred MeV. As in several 
previous works [13, [H-Sl, [MH, [13, [liM [zl this will 
be done in first order perturbation theory. The conse- 
quences of the particle-particle force in the S = 1 , T = 
channel, within the PQRPA will also be examined. The 
importance of this piece of the residual interaction was 
discovered more than 20 years ago by Vogel and Zirn- 
bauer [73 and Cha [7y|, and since then the QRPA be- 
came the most frequently used nuclear structure method 
for evaluating double /3-decay rates. 

A few words will be devoted as well as to the non- 
relativistic formalisms for neutrino-nucleus scattering. 
The most popular one was developed by the Walecka 
group [3J, |77(-[79| , where the nuclear transition matrix el- 
ements are classified as Coulomb, longitudinal, transverse 
electric, and transverse magnetic multipole moments. We 
feel that these denominations might be convenient when 
discussing simultaneously charge-conserving, and charge- 
exchange processes, but seems unnatural when one con- 
siders only the last ones. As a matter of fact, this ter- 
minology is not often used in nuclear /3-decay, ^-capture, 
and charge- exchange reactions where one only speaks of 
vector and axial matrix elements with different degrees of 
forbiddenness: allowed (GT and Fermi), first forbidden, 
second forbidden, etc., types [80, l8l|. There are excep- 



tions, of course, as for instance, is the recent work of Mar- 
ketin et al. [52| on muon capture, where the Walecka's 
classification was used. 

The formahsm worked out by Kuramoto [7J| is also fre- 
quently used for the evaluation of neutrino-nucleus cross- 
sections. It is simpler than that of Walecka, but it does 
not contain relativistic matrix elements, nor is applicable 
for muon capture rates. 

More recently, we have introduced another formal- 
ism [4^-[50] . Besides of being almost as simple as that of 
Kuramoto, it retains relativistic terms and can be used 
for /i-capture. This formalism is briefly sketched here, 
including the consequences of the violation of the Con- 
served Vector Current (CVC) by the Coulomb field. It 
is furthermore simplified by classifying the nuclear ma- 
trix elements in natural and unnatural parities. We also 
show how within the present formalism both the sum rule 
approach, and the formula for ERL look like. 

In Section |ll] we briefly describe the formalism used 
to evaluate different weak interacting processes. Some 
details are delegated to the Appendices: |X| Contri- 
butions of natural and unnatural parity states to the 
transition rates, |B|) Sum rule approach for the inclu- 
sive neutrino-nucleus cross section, [C|) Formula for the 
inclusive neutrino-nucleus cross section at the extreme 
relativistic limit, and|D|) Formula for the muon capture 
rate. In Section lllll we present, and discuss the numeri- 
cal results. Comparisons with experimental data, as well 
as with previous theoretical studies, are done whenever 
possible. Here we firstly sketch the two theoretical frame- 
works, namely the PQRPA and RQRPA, used in the 
numerical calculations. In subsections HIl and IIIIBI we 
present the results for the exclusive and inclusive pro- 
cesses, respectively. Finally, in Section |IV] we give a brief 
summary, and final conclusions. 



II. FORMALISM FOR THE WEAK 
INTERACTING PROCESSES 



The weak Hamiltonian is expressed in the form [771 . [ 
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(2.1) 



relativistic form reads "^ 



Jn = 9v + igA+9pi)cr ■'^ + 9a 
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J = -gAcr - ig^jcr x k - ^^k + gp^icr ■ k)k - g, 
where k = k/|k|. The quantity 



(2.2) 

" M ' 



(2.3) 



is the momentum transfer, M is the nucleon mass, and 
Pi and Pf are momenta of the initial and final nu- 
cleon (nucleus). The effective vector, axial- vector, weak- 
magnetism and pseudoscalar dimensionless coupling con- 
stants are, respectively 



gv — li 5/1 — Ij 9m — '^p Kn — 6.1V, 



gp = gj 



k^ + m^ 



(2.4) 



where the following short notation has been introduced: 



5v 



9v-zr7; .9a = 9A-zrT; ffw = (Sv + gm) 



■2M' 



2M 

9pi = .9p7T77 ; 5p2=5p 



2M' 



2Mmf 



2M mi ' 



(2.5) 



where k = |k|. The above estimates for g^^, and gp come 
from the conserved vector current (CVC) hypothesis, and 
from the partially conserved axial vector current (PCAC) 
hypothesis, respectively. The finite nuclear size (FNS) 
effect is incorporated via the dipole form factor with a 
cutoff A = 850 MeV, i.e., g^g [AV(A2 + k^)] ^ 

In performing the multipole expansion of the nuclear 
operators 



Oa = (O,iO0) =: J„e 



— ik] 



(2.6) 



it is convenient 1) to take the momentum k along the z 
axis, i.e.. 



g^.k.r ^ ^^-V4^(2L+l)jL(p)rLo(f), 
L 

= ^z-V4^(2J + l)jj(p)>0o(f-), (2.7) 
J 

where p = nr, and 2) to define the operators Oq as 

Oa = V4^^rV2J + lOaj. (2.8) 



where G = (3.04545 ± 0.00006) xlO^^^ is the Fermi cou- 
pling constant (in natural units), the Icptonic current 
la = {1, ilfi} is given by the Eq. (2.3) in Ref. [i^l and 
the hadronic current operator J^ = {J,ij0} in its non- 



^ As in Rcf. |49| we use the Walecka's notation |78|| with the Eu- 
clidean metric for the quadrivectors, and o = 1, 2, 3, 4. The only 
difference is that we substitute his indices (0, 3) by our indices 
(0, 0), where we use the index for the temporal component and 
the index for the third spherical component. 



In this way wG avoid the troublesome factor i ■'. In spher- 
ical coordinates (m = —1,0,-1-1) we have 

J$ = .9v + (5a + SpJo'o + ig^M" V • V 

Jrn = -gAf^m + mg^a,n + SmQ[-g^ + ffpaCro] 

- ig^M~^Vra, (2.9) 



and 



Oflj = j^{p)YJo{r)JtD, 
0„u = 5I«^~'^LWL(p)[ii(f)®J]j, (2.10) 



where 



F, 
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(-)J+^V2lTT( 

(-)i+'"(l,-mJm|LO) 



L 1 J 

— m m 



(2.11) 



is a Clebsch-Gordan coefficient. ^ 
Explicitly, from 



O0J = gyM]' +tg^Mf + t{g^+g^,)M-^j, (2.12) 
+ gvM'^j ~ S„,og^M'j . 



(2.13) 



The elementary operators are given by 

M^ = M-'jj{p)YAr){a-V) 



(2.14) 



and 



■M,nJ = E*'"'~'^LWL(p)[il(f)®<T]j, 
L>0 

Ml^j = M"i^zJ-L-ii^uWL(p)[ii(f)®V]j. 



L>0 



(2.15) 



The CVC relates the vector-current pieces of the oper- 
ator dMl) as (see Eqs. (10.45) and (9.7) from Ref. [8l|) 



with 



, = fc0 - 5(A^coui - AM), 



(2.16) 

(2.17) 



* Their values are: 
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2(2J-H)' J.J-i.±i \/2(2J-H)' 



0, i^j,j,±i = T 



1 



where 



Ai^Coul — 



6e^Z 
5R 



= lAbZA-^/^ MeV, (2.18) 



is the Coulomb energy difference between the initial and 
final nuclei, AM = M„ - Mp = 1.29 MeV is the neutron- 
proton mass difference, and S = ±1 for neutrino and 
antineutrino scattering, respectively. 

The consequence of the CVC relation (|2.16p is the sub- 
stitution 

g.M^j - g,M] ^ %.M] , (2.19) 

in (j2.13p . and 0„ij now reads 

+ \m\gyMrnj + SrnO—g^M]'. (2.20) 



The second term in (|2.17p comes from the violation of 
the CVC by the electromagnetic interaction. Although 
it is frequently employed in the nuclear /3-decay, as far 
as we know, it has never been considered before in the 
neutrino-nucleus scattering. AE'coui is equal to 3.8,9.8, 
and 20.0 MeV for ^^C, ^^Fe, and ^ospb, respectively, and 
therefore the just mentioned term could be quite signifi- 
cant, specially for heavy nuclei. 

The transition amplitude for the neutrino-nucleus re- 
action at a fixed value of k, from the initial state |0+) 
in the {Z,N) nucleus to the n-th final state |J^) in the 
nucleus {Z ± 1, iV =F 1), reads 

rj.(^)^5]|(j;j|i/„.(«)|o+)f. (2.21) 

The momentum transfer here is k = pi ~ q^, with pi = 
{pe, iEi\ and q^, = {(\^, iE^}, and after some algebra [4^ 
one gets 



TM^) 



where 
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(2.23) 



are the lepton traces, with 9 = q^-pi being the angle be- 
tween the incident neutrino and ejected lepton momenta, 
and 



go 
Po 



k ■ q,y 
k-pi 



E,{\pe\cose~E,) 

K 

\pe\{\Pi\ - E^cosO) 



(2.24) 



the z-components of the neutrino and Icpton momenta. 
The exclusive cross-section (ECS) for the state |J^), as 
a function of the incident neutrino energy E,j, is 

a({Jl,E,) = MEiF{Z + S,Ee) f d{cos e)Tj^Jn) , 

(2.25) 
where 



Ef 









\pe\ = ^{E^-Ui^J^ 



!-£, 



2E,{Ei 



■ OJ 



J;' 



(2.26) 



and wj^^ = — fcg ^ Ei,~ E^ is the excitation energy of the 
state \]^) relative to the state |0+). Moreover, F{Z + 
S,Ei) is the Fermi function for neutrino {S = 1), and 
antineutrino {S — —1), respectively. 

As it well known the charged-current cross sections 
must be corrected for the distortion of the outgoing Icp- 
ton wave function by the Coulomb field of the daughter 
nucleus. For relatively low neutrino energies (~' 40 — 50 
MeV) this correction was implemented by numerical solu- 
tion of the Dirac equation for an extended nuclear charge 
[8l| . At higher energies, the effect of the Coulomb field 
was described by the effective momentum approximation 
(EMA) [83], in which the lepton momentum pg and en- 
ergy El are modified by the corresponding effective quan- 
tities (see HH, Eq. (34) and (45)]). 

Here, we will also deal with inclusive cross-section 
(ICS), 
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FIG. 1: (Color online) ^^B and ^^N ground state energies (in 
units of MeV) and GT B-values within the PQRPA for differ- 
ent s.p. spaces, as function of the pp-coupling t. The experi- 
mental values are: £;(^^B) = 13.37 MeV, and _B(^^N) = 17.33 
MeV m, and ^(^^B) = 0.466, and ^(^^N) = 0.526 [H. 



ME.) 



J; 



'^ei-^Z^Ei,), 



(2.27) 



as well as with folded cross-sections, both exclusive, 



atiVJ 



dE,aeiJl,E,)ne{E,), 



and inclusive 



CTf = / dE„ae{E^)ni{E„), 



(2.28) 



(2.29) 



where ni(E^) is the neutrino (antineutrino) normalized 
fiux. In the evaluation of both neutrino, and antineutrino 
ICS the summation in ()2.27p goes over all n states with 
spin and parity V < 7* in the PQRPA, and over J^ < 
14± in the RQRPA. 

In the Appendix |A] we show that the real and imag- 
inary parts of the operators OqJ, given by ()2.12p and 
(|2.20p . contribute to natural and unnatural parity states, 
respectively. This greatly simplifies the numerical calcu- 
lations, because one always deals with real operators only. 
In Appendix [D] are also shown the formula for the muon 
capture process within the present formalism. 



III. NUMERICAL RESULTS AND DISCUSSION 

The major part of the numerical calculations have been 
done within the PQRPA by employing the (5-interaction 
(in MeV fm3) 

V^-4TT{vsPs+vtPt)6{r), 

with singlet (vs), and triplet (vt) coupling constants dif- 
ferent for the particle- hole (vh) , particle-particle (pp), 
and pairing {pair) channels [8J]. This interaction leads 
to a good description of single and double /3-decays 
and it has been used extensively in the literature |85l - 
1881 ] ■ The single-particle wave functions were approxi- 
mated with those of the HO with the length parame- 
ter b = 1.67 fm, which corresponds to the oscillator en- 
ergy huj = 45A~i/3 - 25A"2/^ MeV. The s.p. spaces 
82,83, S4, and Sfi will be explored. 

In Rcfs. [4^, |4^, where the 53 space was used, we have 
pointed out that the values of the coupling strengths wf^ 



and vf used in N > Z nuclei (uf 



„PP > 



„pp 



), might not be suitable for N = Z nuclei. In fact, the 
best agreement with data in ■^'^C is obtained for: i) the 
energy of the ground state in ^^N, E{^^N), u) the GT B- 
values in -^^C, B{^^'N) and i?(^^B), and in) the exclusive 
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FIG. 2; Exclusive ^^C{iy,e~y^N cross-section ae{E,,,l^) (in 
units of 10~*^ cm^), plotted as a function of the incident 
neutrino energy _E„. Results for several single-particle spaces 
Sn, and f = 0, within three different energy intervals, are 
shown. The experimental data in the DAR region are from 
Ref. m. 



muon capture in ^^B. A'^'''^ = A'lj'"), is obtained when the 



pp channel is totally switched off, i.e., v^^ 
adopted ph coupling strengths are v^^^ = 



pp 



0. The 
27 MeV fni3 



and v. 



ph 



64 McV fm'^ [481 . For the pp channel it is 



convenient to define the parameters 



,,pp 



patr ■ 



PP 



v: 



pair ' 



where 



(wP"''(p) + «r*''(^))/2 



As in our 



previous work on ^^C, we will use here the same singlet 
and triplet pp couplings, i.e., s = t [43,140]. The states 
with J'^ = 0+, and J'^ = 1+ only depend on s, and t, 
respectively, while all remaining depend on both coupling 
strengths. 

The s.p. energies and pairing strengths for 5*2, 6*3, and 
54 spaces, were varied in a x^ search to account for the 
experimental spectra of odd- mass nuclei ^^C, ^^B, ^'^C, 
and ^^N, as explained in Ref. [431. This method, however, 
is not practical for the space Sq which comprises 21 s.p. 
levels. Therefore in this case the energies were derived in 
the way done in Ref. [5l[ , while the pairing strengths were 
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FIG. 3: (Color online) Same as in Figure [21 but now f = for 
S2, and S3, t = 0.2 for S4, and t = 0.3 for 5*6. The SM, and 
EPT calculations are, respectively, from Refs. [93|, and lol. 
The experimental data in the DAR region are from Ref. 23^ . 



adjusted to reproduce the experimental gaps in ^^C [SOj . 
considering all the quasiparticle energies up to 100 MeV. 

For the purpose of the present study, we also employ 
the RQRPA theoretical framework [90| • In this case the 
ground state is calculated in the Relativistic Hartree- 
Bogoliubov model (RHB) using effective Lagrangians 
with density dependent meson-nucleon couplings and 
DD-ME2 parameterization [91[, and pairing correlations 
are described by the finite range Gogny force 
of the formalism can be found in Refs. [93l[94 



92|. Details 



The RHB 

equations, and respective equations for mesons are usu- 
ally solved by expanding the Dirac spinors and the meson 
fields in a spherical harmonic oscillator basis with ^20 s.p. 
space. In the present study of neutrino-nucleus cross sec- 
tions, with energies of incoming neutrinos up to 600 MeV, 
we extend the number of oscillator shells up to iV = 30 in 
order to accommodate s.p. states at higher energies nec- 
essary for description of cross sections involving higher 
energies of incoming (anti)neutrinos. The number of 2qp 
configurations in the RQRPA is constrained by the max- 
imal excitation energy E2qp- Within the RHB+RQRPA 
framework the oscillator basis is used only in RHB to 
determine the ground state and single-particle spectra. 



The resulting wave functions are converted to coordi- 
nate space for evaluation of the RQRPA matrix ele- 
ments. However, it is the original HO basis employed 
in RHB that determines the maximal £'2qp and the size 
of RQRPA configuration space. 



A. Weak interaction properties of ^^N and ^^B 
ground states 

Let us first compare the QRPA and PQRPA within the 
smallest configuration space 5*2, which contains 16 V = 
f + states, and with null pp coupling: t — 0. The PQRPA 
ground state energies in ^^N, and ^^B, are, respectively: 
a;+i(l+) = 18.319 MeV, and w-i(l+) = 12.528 MeV, 
while the corresponding wave functions read 

^N) = 0.963|lp-/2K/2> 

+ 0.122|lp-/2lP3/2) 



0.2321 lpJ/2lp^/2) 
0.1051 IPJ/2IK/2) 



(3.1) 



and 



^B) = -0.971|lp^/2lp3^/2) 

- 0.125|lpJ/2lK/2> ^ 



+ 0.2041 IPJ/2IP3/2) 
0.0901 IPJ/2IK/2) 



+ ••• (3.2) 

The analogous QRPA energies are quite similar: 
a;+i(l+) = 17.992 MeV, a;_i(l+) = 12.437 MeV. How- 
ever, the wave functions are quite different. The main 
difference is in the fact that QRPA furnishes the same 
wave functions for all four nuclei ^^N, ^°B, ^*N, and ^^B, 
being that of the ground state: 



II 



GS' 



-0.2721 lpJ/2lpr/2) 
0.356|lp^/2lP3/2> - 



-0.759|lpJ/2lp^/2) 
0.4721 lp^/2lK/2> 



+ ■■■. (3.3) 

The difference in the wave functions is an impor- 
tant issue that clearly signalizes towards the need for 
the number projection. In fact, the PQRPA yields 
the correct limits (IP3/2 -^ ^Pi/2 ^^'^ ^P3/2 ~^ ^Pi/2) 
for onc-particle-one-hole (Iplh) excitations on the ^^C 
ground state to reach the -'^^N, and ^^B nuclei. All re- 
maining configurations in p.ip . and p.2p come from 
the higher order 2p2h, and 3p3h excitations. Con- 
trary, the QRPA state (|3.3p is dominantly the two- 
hole excitation [{lp'^,^)~^, {lp'^,^)~^], which corresponds 

to the ground state of ^°B. More details on this ques- 
tion can be found in Figure 3 of Ref. ^^. The Iplh 

amplitudes[(lp3,2)^^, 1^1/2]' ^^'^ [(1^3/2)"^' (^^1/2)] ^^'^ 
dominantly present in the following QRPA states 



II2") = 



ii,r) = 



0.708|lp^2lP3/2>+0-703|lpJ/2lK 



3/2 



1/2/ 



-0.4761 lpJ/2lK/2> + 0.4371 lpJ/2lp^/2) 
0.441|lp-/2lP3/2> - 0.0961 K/2lpr/2) 
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^Cii>,e+y'^'B cross- 
"^^ plotted as a func- 



(3.4) 



FIG. 4: (Color online) The calculated 

section a^+ {Ep, 1^) (in units of 10~*^ cm" 

tion of the incident antineutrino energy _Ef . Same as in Figure 

O the value of i is for the s.p. spaces 52, and S3, 0.2 for 

54 and 0.3 for Se- The EFT calculation from Ref. [Q is also 

shown. 

The wave functions displayed above clearly evidence 
the superiority of the PQRPA on the QRPA. Therefore 
from now on only the PQRPA results will be discussed 
for the exclusive observables. 

In Figure [1] we show the ^^B and ^^N ground state 
energies, and the corresponding GT i?-values within the 
PQRPA for different s.p. spaces, as function of the pp- 
coupling t. One sees that the energies depend rather 
weakly on both, and agree fairly well with the measured 
energies: E{^^B) = 13.37 MeV, and E^-^N) = 17.33 
MeV [931, although the first one is somewhat under- 
estimated, while the second one is somewhat overesti- 
mated. Both GT S-valucs significantly increase with t 
and diminish when size of the s.p. space is increased. For 
spaces 5*2 and S3 the best overall agreement with data 
(S(12b) = 0.466, and ^(i^N) = 0.526 ^) is achieved 
with t ~ 0, while for spaces ^4 and Sq this happens when 
the couplings are, respectively, t = 0.2, and t = 0.3. 

After establishing the PQRPA parametrization, we an- 
alyze the behavior of the ECSs of the ground states in 
^^N and ^^B, as a function of the size of the configura- 
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FIG. 5: (Color online) Muon capture transition rate to the 
^^B ground state in units of 10^s~^, and electron and muon 
folded ECSs to the ^^N ground state in units of 10~*^ cm^ and 
10^^^ cm^, respectively. The experimental values, in above 
units^re: A(^^B) = 6.2 ±0.3 ^, aeC^) = 9.1 ± 0.4 ± 
0.9 [23, and aeJ^N) = 8.9 ± 0.3 ± 0.9 H, and auC^N) = 
6.6 ± 1.0 ± 1.0 [ill, and ^^(^^N) = 5.6 ± 0.8 ± 1.0 [H]. 



tion space. Figure [2] shows the ECSs for the reaction 
^^C{i',e~y^N (in units of 10^^^ cni^) for several config- 
uration spaces, and for t = 0, within three different en- 
ergy intervals. The top panel represents the DAR region, 
where experimental data are available [23 , and search for 
neutrino oscillations was done [25J, 127[ . The middle panel 
represents the region of interest for supernovae neutrinos, 
as pointed out in Refs. [la, [931' while the bottom panel 
shows the asymptotic behavior of the cross-section, which 
becomes almost constant for E^ ~ 200 MeV. Within the 
spaces ^2 and ^3 the calculations reproduce quite well 
the experimental cross sections in the DAR region, as 
seen from the first panel. 

In Figure [3] we show the calculated ECSs for the re- 
action ^^C(i', e^)^^N within several configuration spaces, 
but now with different values of the pp-coupling. From 
comparison with the experimental data in the DAR re- 
gion [23 one observes that the appropriate values for the 
coupling t for s.p. spaces ^4, and Sq, are, respectively, 
t = 0.2, and t = 0.3, i.e., the same as those required to re- 
produce the experimental energies and the GR _B-values 
in i^B, and ^^N. 
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FIG. 6; (Color online) Inclusive ^^C{u,e~y^N cross-section 
(7g-(£'^)(in units of 10"'^^ cm^) plotted as a function of the 
incident neutrino energy E,,. The PQRPA results within the 
s.p. spaces S2, S3, and 56, and the same values of s = i 
as in Figure |31 are compared with two sum rule limits (as 
explained in the text): SRgt, and SRif obtained with aver- 
age excitation energy cJjT of 17.34, and 42 MeV, respectively. 
Several pre vious RPA-like calculations, namely: RPA ij. 
CRPA [io2, and RQRPA within S20 for S29p=100 MeV [5l|, 
as well as the SM [431, and the TDA [3J| are also shown. 



This change of parametrization hint at the self- 
consistency of the PQRPA, and comes from the fact 
that in this model: i) the GT strength allocated in the 
ground state is moved to another 1+ states when the 
size of the space is increased, and ii) the effect of the 
pp residual interaction goes in the oppositive direction, 
returning the GT strength to the l^*" state. Only for 
the space 52 the cross-section a^- {E^, 1]^) is appreciable 
larger (at i?^ > 60 MeV) than for other spaces, which 
is just because of the small number of configurations in 
this case. In the same figure are exhibited as well the 
results for the ECSs evaluated within the SM [9^, and 
the EPT [1^ . Both of them agree well with the data and 
with the present calculation. 

The results for the reaction {i>,e^) to the ground 
state in ^^B are shown in Figure SI The cross-section 
cre+ {Ej}, IJ") is similar to that produced by neutrinos but 
significantly smaller in magnitude. When compared with 



10 



the EPT result [l6[, which are also shown in the same 
figure, one notices that they are considerable different. 
To some extent this is surprising as in the case of neutri- 
nos the two models yield very similar results. One should 
remember that in the EPT model the axial form factor, 
used for both neutrinos and antineutrinos, is gauged to 
the average of the GR i?-values in ^^B, and ^^N, which, 
in turn, arc well reproduced by the PQRPA. Therefore 
it is difficult to understand why the EPT results agree 
with the present calculations for neutrinos and disagree 
for antineutrinos. 

In Figure [5] we show the dependence on the con- 
figuration space of the exclusive muon capture transi- 
tion rate A(l^) to the ^^B ground state, and the elec- 
tron and muon flux-averaged ECSs, given by (|2.28p . 
to the -^^N ground state, i.e., a'e(lf), and a^{li). 
As in Refs. |4^ |4^ the electron neutrino distribution 
rie (-El/) was approximated with the Michel energy spec- 
trum |9l.l9SJ|. and for the muon neutrinos we used nfj,{E^) 
from Ref. [30|- The energy integration is carried out in 
the DAR interval rue+coj^ < A^^^ < 52.8 MeV for elec- 
trons and in the DIE interval m^ 4- ujjf < Aj^^ < 300 
MeV for muons. Erom Eigure EJ and comparison with 
experimental data: 

A(12b) = 6.2 ±0.3 m, 

CTe(i2N) = 9.1 ± 0.4 ± 0.9 P[, 8.9 ± 0.3 ± 0.9 M, 

ct^(12n) = 6.6 ± 1.0 ± 1.0 lH, 5.6 ± 0.8 ± 1.0 [H, 
one finds out, as for results shown in Eigurcs[Tl and|3l the 
model self-consistency between s.p. spaces and the pp- 
couplings. That is, for larger s.p. spaces larger values of t 
are required. In brief, the experimental data of CTe(^^N), 
and CT^(^^N) are well reproduced by the PQRPA. The 
same is true for the SM calculations [42, |43], while in 
RPA, and QRPA models they are strongly overestimated, 
as can be seen from Table II in Ref. [43, and Table 1 in 
Ref. [si. 



B. Inclusive cross-sections ^^C{u,e )^^N and 

^^C{u,e+y^B, and Sum Rule 



In Eigureiniwe confront the PQRPA results for the ICS 
a^-{E^) within spaces S2, S^, and Sq with the corre- 
sponding sum-rules (jf3^{E^) evaluated from (|Bip . One 
immediately sees that the PQRPA results depend very 
strongly on the size of the employed s.p. space. On the 
other hand, as already mentioned in the Appendix |B1 the 
sum rule crf?'{E^) depends on the average energy wj^. 
Here we use two values ZUj^ = 17.34 MeV, which is the 
ground state energy ^^N, (GT-rcsonance), and aJjT = 42 
MeV, whic h is roughly the energy of the first forbidden 
resonance |100l | . The corresponding curves in Eigure [5] 
are labeled, respectively, as SRgt, and SRif- They 
should be the upper limits for allowed and first forbidden 
transitions, respectively. The validity of these sum rules 
is questionable for neutrinos energies of several hundred 
MeV, as already pointed out by Kuramoto et al. [2J| . In 
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FIG. 7; (Color online) Inclusive ^^C(t', e^)^^B cross-section 
a^+{E„) (in units of 10~^® cm^) plotted as a function of the 
incident neutrino energy Ea. All results were obtained in the 
same way as in the neutrino case in Figure (6] 



fact, we note that the cross section SRgt {SRif) exceeds 
the free particle cross section erg = &o_(jye -t- n — > e~ + p) 
for E^ > 200 MeV {E^ > 300 MeV) [lOlj . 

Several previous SM and RPA-like calculations of 
tTe-(Si/), employing different effective axial-vector cou- 
pling constants, and different s.p. spaces, are exhibited 
in Figure El as well, namely: 

a) TDA [SJ, witlio4 = 1.23, and S2, 

b) SM and RPA ^, with g^ = 0.88, and S3, 

c) CRPA [lOl, with g^ = 1.26, and S4, 

d) RQRPAJHll, with g^ = 1.23, S20, and £;2qp=100 MeV. 
It is important to specify the values 5^ because the par- 
tial cross sections arc predominantly of the axial- vector 
type (specially the allowed ones) , which are proportional 
to g^ . In spite of very significant differences in g^ , and 
the s.p. spaces, the different calculations of o-g- (E^) yield 
quite similar results for energies E^, < 130 MeV, lying ba- 
sically in vicinity of the the sum-rule result SRif- But 
for higher energies they could become quite different, and 
arc clearly separated in two groups at Ei, ~ 300 McV. In 
the first group with (T^-{E,y) ^ 5 x 10"'^^ cm^ arc: the 
SM, TDA, and PQRPA within spaces S2, S3, while in the 
second one with a^- {Ey) > 10 x 10^^^^ cm^ are: the RPA, 
RQRPA, CRPA and PQRPA within space S'e. Volpe et 
al. [43I have found that the difference between their SM 
and RPA calculations is due to differences in the corre- 
lations taken into account, and to a too small SM space. 
We also note that only the CRPA result approaches the 
sum rule limits for E^, > 200 MeV. 

Similar results for the inclusive ^'^C{v,e^Y^'Q cross- 
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FIG. 8: (Color online) Inclusive ^"^Giy^ e~)^'^N cross-section a^- (_E„)(in units of 10^^® cm'^) plotted as a function of the incident 
neutrino energy Eu^ evaluated in RQRPA with different configuration spaces. These cross sections are plotted as functions of 
the incident neutrino energy with different cut-off of the E2q^ quasiparticle energy as it is explained in the text. The left and 
right panels show the cross section evaluated with S20, and S30 s.p. spaces. The last cross section shows that the convergence 
of the calculation is achieved up to 600 MeV of incident neutrino energy. 



section a^+(Ej^) arc displayed in Figure [TJ and analo- 
gous comments can be done here. For the comparison, 
wc show in the figure the antineutrino-^^C cross-sections 
evaluated with the CRPA }102j . 



C. Large configuration spaces 

As there are no experimental data on flux unfolded 
ICSs for E,y < 400 MeV we cannot conclude which of the 
results displayed in Figures [6l and [7] are good and which 
are not. We can only conclude that the ICSs strongly 
depend on the size of the s.p. space. In the PQRPA 
calculations wc are not able to use spaces lager than Sq 
because of numerical difficulties. Thus instead of using 
the PQRPA, from now on we employ the RQRPA where 
such calculations are feasible. It is important to note 
that within the RHB-I-RQRPA model the oscillator basis 
is used only in the RHB calculation in order to determine 
the ground state and the single-particle spectra. The 
wave functions employed in RPA equations are obtained 
by converting the original HO basis to the coordinate 
representation. Therefore, the size of the RQRPA con- 
figuration space and 2qp energy cut-offs are determined 
by the number of oscillator shells in the RHB model. 

First, we analyze the effect of the cut-off energy within 
the 5*20 space on a^- {E„) for E^ up to 600 MeV. From 



the left panel in Figure [S] one sees that at high energies 
this cross-section increases roughly by a factor of two 
when E2qp is augmented from 100 to 200 MeV. The in- 
crease of the cross-section is also quite important when 
E2qp is moved from 200 to 300 MeV. For the limiting 
value of i?2qp=300 MeV, all possible configurations are 
included in RQRPA calculations. Next, we do the same 
within the 5*30 space, and the resulting (j^-{Ey) are dis- 
played on the right panel of Figure HI From the com- 
parison of both panels it is easy to figure out that up 
300 MeV the cross sections obtained with the 



to E' 



2qp 

S30 space are basically the same to those calculated with 
the 5*20 space. Small differences between the cross sec- 
tions using S20 and 6*30 spaces for E2qp up to 300 MeV 
are caused by modifications of positive-energy single- 
particle states contributing to the QRPA configuration 
space within the restricted 2qp energy window. But, for 
El, > 400 MeV additional transition strength appears in 
the 5*30 space when E2qp is moved up to 400 MeV, from 
where further increase in E2qp has a very small effect. 
We conclude therefore that the configuration space en- 
gendered by TV = 20 HO shells with E2qp = 300 MeV, 
is large enough to describe cr^-{Ei,) with E^ up to 400 
MeV. Similarly, the space brought about by iV = 30 HO 
shells with E2qp = 400 MeV is appropriate to account 
for cr^-{E^) up to Ei, — 600 MeV. For larger neutrino 
energies very likely we would have to continue increasing 
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FIG. 9: (Color online) Same as in Figure[8]but for the ^^C{u,e'^)^'^'B cross-section a^+{Ec,). 



the number of shells. 

Analogous results for antineutrino ICSs ae+{Ep) are 
displayed in Figure [9l One notes important differences 
in comparison with ae-{E,j) shown in Figure [51 First, 
here the spaces 5'2o and 6*30 yield almost identical results 
in the entire interval of antineutrino energies up to En ~ 
600 MeV. Second, the successive increase in the cross- 
sections when the cut-off E2qp is augmented in steps of 
100 MeV are smaller, and decrease more rapidly than in 
the neutrino case. This suggests that the configuration 
space is now sufficiently large to appropriately account 
for (jg^+ (Eii ) even at antineutrino energies larger than 600 
MeV. 

At present, due to numerical difficulties, we cannot 
perform the RQRPA calculations for the full range of 
neutrino energies where the QE cross section was mea- 
sured at MiniBooNE [13], but only up to 0.6 GeV. How- 
ever, we feel that it could still be illustrative for com- 
parison with data. This is done in Figure IIOI which 
is basically a piece of Figu re 21 Ref. [23 for the QE 
a^- {Ey) (see also Ref. [l03j ). where is incorporated our 
result for a^- [E^) from Figure [8] for 5*30 and E2qp = 500 
MeV. As already pointed out in the Introduction, at rel- 
atively high energies {Ey > 300 MeV) the electron and 
muon neutrino cross sections converge to each other, and 
therefore, in the present analysis, the electron neutrino 
cross section provides a reasonable upper limit estimate. 
One sees that we underestimate the data by almost a 
factor of two. But one should keep in mind that, while 
we use Qa = 1 (see (I2.4p) in the RFGM calculation done 
by Martini et al. [73, Il03| Qa = 1.255 was used. Be- 



ing the axial-vector contribution dominant for the latter 
value of the coupling constant, one would have to re- 
normalize our Uf,- {Ey) by a factor ~ 1.5. Such a result 
is also shown in Figure [TUl and although the resulting 
cross section still underestimates somewhat the data for 
a^-{Eu), it is notably superior to the pure Ip-lh result 
from Ref. [Tj], where good agreement with the data is 
achieved only after considering additional two-body (2p- 
2h) and three-body (3p3h) decay channels. One should 
keep in mind, however, that as the weak decay Hamilto- 
nian is one-body operator, these excitations are only fea- 
sible via the ground-state correlations (GSC), which ba- 
sically redistribute the Ip-lh transition strength without 
increasing its total magnitude when the initial wave func- 
tion is properly normalized. In the present work, as well 
as in all SM-like calculations, the GSC, and a normalized 
initial state wave function are certainly considered to all 
orders in perturbation theory through the full diagonal- 
ization of the hamiltonian matrix. On the other hand, in 
Refs. |73l . Il03l | the GSC are taken into account in second 
order perturbation theory, but there are no references to 
the normalization of the ^■^C ground-state wave function. 
How to carry out the normalization in the framework of 
the infinite nuclear matter model is discussed in a re- 
cent paper related to t he nonnicsonic weak dec ay of the 
hypernucleus ]^C [lo3; see also Refs. [loMo3. 
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FIG. 10: (Color online) The calculated RQRPA (within S30 
and E2qp — 500 MeV) quasi-elastic (i/g/'^C) cross section per 
neutron (full line) is compared with that for the {Ufij'^^C) scat- 
tering data measured at MiniBooNE [IJ; with dotted line 
is shown the same calculation but renormalized by a factor 
1.5. Al so a re displayed the calculations done by Martini et 
al. [z3, [10^ within the RFGM for pure (Ip-lh) excitations 
(dashed line), and with the inclusion of the np-nh channels 
(dot-dashed line). 



D. Multipole decomposition of cross-sections 

We did not mention yet the contributions of different 
multipoles to the ICSs. Normally, the RHB model within 
520, and with J'^ < 7 , provides converged results for 
RQRPA excitation spectra at incident neutrino energies 
E^, < 300 MeV as seen from Figure 2 in Rcf. [Hj. But 
this is not the case for neutrino-nucleus cross sections at 
energies Ei, > 300 MeV where one has to consider large 
cutoff energies i?2gp- In fact, it is necessary to consider 
more and more multipoles according as the configuration 
space is enlarged by increasing E2qp. This is illustrated 
in Figure [TT] for the case of E2qp = 500 MeV. One sees 
that are significant all multipoles up to J'^ = 14* for 
neutrinos, and up to J'^ = 11 for antineutrinos. 

Next we discuss the partial multipole contributions to 
the ICS, having in view the degree of forbiddenness of 
the transition matrix elements, namely, 

• Allowed: cr^+{Eo), with V = 0+, 1+, 

• First-forbidden crlf^iEp), with V = Q- , 1^,2", 

• Second-forbidden crl^{E^), with J'^ = 2+,3+, and 

• Third-forbidden crff (S^) with J'^ = 3", 4", 
cross-sections. Thus, in the left panel of Figure [12] 
we show these individual contributions for the inclusive 



C(j>, e"'")^^B cross-section ag+{Ep), evaluated within 
both the PQRPA (spaces ^2, and Se) the RQRPA (space 
530 with E2gp = 500 MeV). 

The same is done for the corresponding derivatives, i.e., 
the spectral functions di7e+ {Eo)/ dEo, on the right panel 
of the same figure. Several conclusions can be drawn. 
First, as in the case of total ae+{Ec,), they depend very 
strongly on the size of the configuration space. This de- 
pendence, in turn, increases with the degree of forbid- 
denness; that is, it is more pronounced for first-forbidden 
than for allowed transitions, and so on. Second, within 
the PQRPA the allowed cross-section af+ {E^) exhibits a 
resonant pattern at low energy, and is dominant for Ej) < 
50 MeV. For large s.p. spaces its contribution is quite sig- 
nificant even at Sp = 500 MeV. ^ In the case of RQRPA, 
the spectral function d(Tf+[Ei))/dEi) also displays low- 
energy resonant structure, and a\ {Ey) is always smaller 
in magnitude than in the PQRPA case. Third, (t^+ (£'£>) 
is peaked at Ec, ^ 75 MeV, and its contribution is always 
larger than that of (T^+{Ei:,) for Eo > 150 MeV. Fourth, 
a'^+{E;}), and a^+{Eii) mainly contribute in the interval 
150 < E'ii < 400 MeV, and their overall contributions are 
of the same order of magnitude, and comparable to that 
of the al+{Ec,). Fifth, the contributions of the remain- 
ing multipoles with J'^ = 4+,5*,6*,7* are always very 
small for the space 6*2, but are quite sizeable for Sq at 
high energies. For instance, at E^ = 100, 300, 600 MeV 
they contribute, respectively with 0.02%, 0.86%, 1.18% 
within the space S2, and 0.04%, 14%, 20% for 5*6. With 
further increase of the single-particle basis, configura- 
tions from higher multipoles become more pronounced 
at higher neutrino energies. In particular, the sum of 
contributions coming from ^'^ =4+,--- ,11* multipoles 
when evaluated within the RQRPA using the space 6*30 
and maximal value of i?2qp=500 MeV, are 1.1%, 14.4%, 
and 33.2% at Ey=100, 300, and 600 MeV, respectively. 

Recently Lazauskas and Volpe |l08L Il09l | have sug- 
gested the convenience of performing nuclear struc- 
ture studies using low energy neutrino and antineutrino 
beams. Because of feasibility reasons the flux covers 
80 MeV only. Nevertheless, from the analysis of ^^O, 
56pg^ iooMo and ^ospb nuclei within the QRPA using 
the Skyrmc force they were able to disentangle the mul- 
tipole distributions of forbidden cross-sections, showing 
that the forbidden multipole contribution is different for 
various nuclei. In this work we extend this kind of study 
to 12c. 

In Table U we show the results for the flux-averaged 
cross sections ag+ for the reaction ^^C{C',e'^y^'B. In 
(|2.29p we have used the same antineutrino fluxes 71^+ {Ej}) 



^ The denominations here don't have exactly the same meaning as 
in the low-energy /3-dccay, where allowed transitions are those 
within the same HO shell (AA^ = 0), while here all values of AA'^ 
are permitted. Similarly happens with the forbidden transitions. 
The degrees of hindrance basically come from value of the orbital 
angular momenta. 
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FIG. 11: (Color online) Left and right panels show, respectively, the cross sections a^-{E^), and a^+{Ep) (in units of 10^" 
cm^) evaluated in RQRPA for S20, and Sgo s.p. spaces with the cutoff -£2,3^ ~ 500 MeV, and different maximal values of J 
with J going from 1 up to 14 for neutrinos, and from 1 up to 11 for antineutrinos. 



as in Rcf. |108j . i.e., the DAR flux, and those produced 
by boosted ^Hc ions with different values of time dilation 
factor 7 = l/^/l — v^ jc?. Results of two calculations are 
presented: i) PQRPA within Se, and ii) RQRPA within 
A^ = 20, and cutoff E-i^j, = 300 McV. One sees that in 
both models, and principally in the PQRPA, the allowed 
transitions dominate the forbidden one, and specially for 
the low-energy beam with 7 = 6. The contributions of 
the second-forbidden processes are very small in all the 
cases, while those coming from third- forbidden ones are 
always negligible. All this is totally consistent with the 
results shown in Figure [T^ from where it is clear that to 
study second and third forbidden reactions in ^^C, one 
would need fluxes 71^+ {En) with Eo at least up to > 150 
MeV. It should also be stressed that our results both for 
allowed and forbi dden transitions fully agree with those 
obtained in Ref. jl08l |: the difference in ^^O comes from 
the double-shell closure in this nucleus. 



E. Supernova neutrinos 



TABLE I: Fraction (in %) of flux-averaged cross sections 
a^+ for ^^C{u,e+y^B for allowed (A), first-forbidden (IF), 
second-forbidden (2F), and third-forbidden (3F) proc esses . 
The antineutrino fluxes n^+ {Ec) are the same as in Ref. [lOSi ]. 
i.e., the DAR flux, and those produced by boosted ®He ions 
with different values of 7 = 1/a/1 — w^/c^. Results of two cal- 
culations are presented: i) PQRPA within 55, and ii) RQRPA 



within A'^ = 30, and cutoff E2q 



■ 300 MeV. 
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We also address briefly the v/v-^'^C nucleus cross- 
sections related with astrophysical applications, the 
knowledge of which could have important implications. 
For this purpose, are evaluated the a^± folded with 
supernovae v/v spectra represented by a normalized 
Fermi-Dirac distribution with temperatures in the inter- 
val T^^ =2—12 MeV, which includes the most commonly 



used values are T^, = 3.2 MeV and Tp^ = 5.0 MeV. For 
mean e nerg i es jE ^) « 3.15 x T^, and zero chemical po- 
tential jllOl llll| the neutrino flux is 



ne{Ey) 



0.5546 Ei 



(3.5) 
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FIG. 12: (Color online) Left panel: Allowed (J" = 0^ 
and third-forbidden (J'^ = 3~,4~) inclusive ^^C{i',e~^ 



, 1+), first-forbidden {J" = 0", 1",2"), second-forbidden {J" = 2+,3+) 
I^^B cross-section a^+ (Ep) (in units of 10 



*^ cm^), plotted as a function 



of the incident neutrino energy Ei^. Right panel: Same as left panel but now for da^+{E[,)/dEo (in units of 10 cm MeV 



and similarly for antineutrinos. For the sake of simplicity 
we do not analyze same relevant aspects of ne{E^) in 
supernova sim ulatio n, such as the MSW effect (see, for 
example. Rcf. |ll2l | ), an d the spectral swapping of the 
neutrino flux (Ref. J113l |). In Figure [13] we confront the 
v—^^C cross sections averaged over supernova I'-fluxes for 
the range of T^ = 2 — 12 MeV, obtained within following 
calculations: 



i) PQRPA within S'e, 

ii) RQRPA within Sy,Q and E2qp = 500 MeV, and 

iii) SM done by Suzuki et al. ^ with the SFO Hamil- 
tonian (the PSDMK2 interaction yields a quite similar 
result). 

As seen from Figure 1131 in the vicinity of the temper- 
atures mentioned at the beginning [T^ = 3 — 5 MeV), 
these three calculations yield, respectively, that: i) CTe+ 
is significantly larger than a^- , ii) a^- is only slightly 
larger than a^+ , and iii) CTe+ = a^^ . Both SM cross sec- 
tions are always smaller than those obtained in the the 
other two calculations, and specially in comparison with 
the RQRPA one. 



IV. SUMMARY AND CONCLUDING 
REMARKS 



The present work is a continuation of our previous 
works [i^li^. In fact, the formalism for weak interaction 
processes introduced there is now elaborated more thor- 
oughly yielding very simple expressions for the transition 
rates, which greatly facilitate the numerical calculation. 
This is done through the separation of the nuclear matrix 
elements into their real and imaginary parts, which, in 
turn, permits to split the transition rates, for neutrino- 
nucleus reactions (|2.22p into natural (jA6p and unnatural 
parity (jA7p operators. Similar separation is done for the 
muon capture transition rate (jD3|) in Eqs. (jD4[) and (jD5p . 
Moreover, consequences of explicit violation of CVC hy- 
pothesis by the Coulomb field (|2.18p are addressed for 
the first time, and the sum rule approach for the in- 
clusive cross section, proper to the present formalism, 
has been worked out in the Appendix B. For the sake of 
completeness, the extreme rclativistic limit of neutrino- 
nucleus cross section is also presented in the Appendix 
C, where in the formula for transition rates turn out to 
be still simpler. We note that, except at very low neu- 
trino energies, they can be used without any restriction 
in practical applications. 

We have discussed in details the inclusive properties 
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FIG. 13: (Color online) Flux-averaged neutrino and antineu- 
trino cross sections a^± in ^ C with typical supernovae fluxes. 



that comprise: 

i) Ground state energies in "'^^B and "'^^N, and the cor- 
responding GT B-values (Figure [Ij, 

ii) Exclusive ^^C(i^, e")"'^^N cross-section ae{Ei,, 1^), as 
a function of the incident neutrino energy E^, (Figures [2l 
and [3]), 

iii) Exclusive ^^C(i>,e"'')^^B cross-section ae+{Ep, if), 
as a function of the incident antineutrino energy Ep (Fig- 
ure m, and 

iv) Muon capture transition rate to the ^^B ground 
state, and electron and muon folded cross-sections to the 
i^N ground state a^ilt), and o'f^ilt) (Figure E]). 

Special attention was paid to the interplay between the 
size of the configuration space, and the magnitude of the 
residual interaction within the pp-channel. It was found 
that as the first becomes larger, the second has to increase 
to obtain the agreement with the experimental data for 
the exclusive observables. 

The main purpose of our discussion of exclusive prop- 
erties was to put in evidence the limitations of the RPA 
and the QRPA models. The basic problem in the imple- 
mentation of the RPA is the lack of pairing correlations, 
i.e. the inability for opening the IP3/2 shell, while defi- 
ciency of the standard QRPA is in the non-conservation 
of the number of particles, as evidenced by the wave func- 
tions ((XT|) . (133), and (1331) presented in Sec. IIIl In this 
way we have definitively established that the SM and the 
PQRPA are proper theoretical frameworks to describe 



the ground state properties of ^^B and ^^N. ^ 

The inclusive cross-sections ^^C(i^, e~)^^N and 
i2C(£>,e+)i2B have been studied within the PQRPA 
in the same manner as the exclusive ones for Ei,^ up 
to 300 MeV. As there are no experimental available 
in this case the comparison is done with the previous 
calculations only ^, and displayed in Figures IHl and 13 
Here, unlike within Figures lU l3j and IH we also show 
the results obta ined with the other RPA-like models 
[SJ, S, ImI, Il02j . which could be a suitable framework 
for describing global nuclear properties such as it is the 
inclusive cross-sections. 

When the size of the configuration space is enlarged 
the calculated PQRPA cross-sections, at difference with 
the exclusive ones, steadily increase, and particulary for 
neutrino energies larger than 200 MeV, in spite of includ- 
ing the particle-particle interaction. At low energy they 
approach to the cross section of the first-forbidden sum- 
rule limit, but are significantly smaller at high energies 
both for neutrino and antineutrino. 

The largest space that we can deal within the number 
projection procedure is the one that includes all the or- 
bitals until the A^ = 6 HO shell. This is the reason why 
we have recurred to the RQRPA where it is possible to 
employ larger configuration spaces. It seems that when 
the number of shells is increased to A^ = 30, and the cut- 
off energy i?2gp is large enough, the cross sections very 
likely converge as shown in Figures lU HI HOl and ITTJ 

The Figure HUl also indicates that the RQRPA is a 
promising nuclear model to reproduce the quasi elastic 
(i'^,^^C) cross section in the region of £'1, ~ 1 GeV which 
has been measured recently at MiniBooNE [1J|. We do 
not know whether the discrepancy between the experi- 
ment and the theory comes from the non completeness 
of the configuration space or from the smallness of the 
effective axial-vector coupling constant that we are using 
Pa = 1- It could also happens that we need ^a = 1 for 
the low energy exclusive cross section and ^a = 1-255 for 
the high energy inclusive cross section. We do not under- 
stand the reason for such a energy depen denc e of g^, but 



it is consistent with the Eq. (23) in Ref. Ill4l | where it is 
shown that for the low energy /3-decay ^a could be much 
more quenched that the total GT strength. We hope to 
be able to say more on this matter in the next future. 

We have also addressed the issue of multipole com- 
position of the inclusive cross sections, by separating 
them into allowed (J'^ = 0"*",1+), first-forbidden (J'^ = 
0-,l-,2-), second-forbidden {J'" = 2+,3+), and third- 



After our work has boon finished, Cheoun et al. 111511 have pre- 
sented a new evaluation of the ECS in ^-^C within the QRPA. 
They get good agreement with data for CTe(^^N), which is at 
variance with the previous QRPA calculation 1431 . 

As already mentioned in the introduction the only available 
experimental data on ^■^C ICS is the low neutrino energy (Ei/^ < 
60 MeV) folded one, which has already been discussed in our 
previous works [49l45lll . 
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forbidden {J^ = 3,4) processes. The results for the 
antiiieutrino reaction ^^C(i>, e"'")^^B are displayed in Fig- 
ure [12] both for the PQRPA and the RQRPA. Of course, 
similar results are obtained also for neutrinos. We re- 
mark that the spectral functions da'^j^{Ei^)/ dE^, when 
evaluated within the PQRPA, clearly put into evidence 
the resonant structure of the allowed cross-section, which 
is mainly of the GT type. 

The study of the partial ICS s has been related with 
the proposal done in Refs. |l08l Il09| of performing nu- 
clear structure studies of forbidden processes by using 
low energy neutrino and antineutrino beams. From the 
results shown in Table U for the flux-averaged cross sec- 
tions CTe+ in the reaction ^^C(i>, e"'")^^B we show that the 
contribution of allowed transitions decreases gradually in 
favor of the first forbidden transitions according with the 
increase of 7-boost. We conclude that to study high for- 
bidden reactions one would need f/-fluxes with E^ up to 
> 150 MeV in ^^C. 

At the end we considered possible astrophysical ap- 
plications of the v/v-^'^C nucleus folded cross sections 
CTe± , using supernovae v/v spectra represented by a 
normalized Fermi-Dirac distribution with mean energy 
{E,j) « 3.15 X T;^, and zero chemical potential. It is found 
that for temperature T^, = 3 - 5 MeV both the PQRPA 
and RQRPA models yield significantly larger cross sec- 
tions that the previously used shell model. 
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are real, A^±ij and A^±ij arc not. Explicitly, 






(A2) 



where 



Mtf ^ M-lf, 



= E i-)^'-^-'^^' Fuudp) [Yl{v) ® <t]j , 



L=J±1 

tA.I 






M''Jl, = -Fw]^{p)[Y^{v)®(T]y 



'ij 



ij 



= E i-)'^'-^-'^'^ Fuijiif) [iL(f) ® V]j 



L-J±l 

1V.1 



Mi; ^ -Air-{j = -^ijjjj(p)[rj(f)®V]j, (A3) 
with L > 0, and J 7^ 0. Thus 

0±ij = t{-g^±g„){Mtf±tMtj') 

+ .gv(A^ir±*A^r/), (A4) 

and writing 

OflJ = Ow +*O0j> 

OrnJ = O^j+lO'^j, (A5) 

it is not difficult to discover that: 

• For natural parity states , with n = (— ) , i.e., J^ = 
0+,l-,2+,3-,---: 

01 = 9vM], 

0|ij = i±g.-g^)Mt,'+gvM';f, (A6) 

• For unnatural parity states, with tt = {—)'^^^, i.e., 
J-" = 0-,l+,2-,3+,---: 

Ogj = -9AM^-{g^+g^,)M^i, 

Ooj = (ffA-flp2)-M0J' 

Oiij = (5,T5w)-Mu"TgvA^r/- (A7) 



Appendix A: Contributions to 7j^ («;) of natural and 
unnatural parity states 

The real and imaginary parts of the operators Oqj 
given by p.l2[) and p.20p do not contribute simultane- 
ously. In fact, the jKGqj (30qj) contributes to natural 
(unnatural) parity states, which means that we always 
can work only with real operators, which greatly simpli- 
fies the calculations. To see this we note that, while the 
operators M]' , Mj, and 



L=J±1 



(Al) 



These operators have to be used in (|2.22p . instead of 
those defined in ([2T2)) . and ([2:20)) . 

The correspondence between the individual matrix el- 
ements, defined by Donnelly, and Peccci in Eq (3.31)of 
Ref. [731, ^''^d the ones used here, is: 

Mj ^ M], 
Aj ^ V2M'^j\ 

Ej ^ V2Mfj\ (AS) 

S^ ^ -V2Mff, 
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Moreover, the correspondence between the Hnear com- 
binations of these matrix elements defined in |77l . Eqs. 
(3.32)] (for Lj see [zi, Eq. (14)]), and the ones intro- 
duced here is: 

• For natural parity states : 



Mi = 0(,j, 
L} — Ooj, 

• For unnatural parity states: 

m! = Ow, 

-it^j = Ooj, 
i(ff5±f7"S) = V20±ij. 

The following relation can also be useful: 
Oflj = TWj, 

OrnJ = ' 



(A9) 



(AlO) 



(All) 



Cj, 



for TO = 



1 

\/2 



mTj^ + Tf 



for TTt = ±1 



where Mj = Mj + M°, Cj = Lj + L], 7]°' = Tf + Tf^, 



mag5 



and Tj'' = Tj""^** + Tj . 

The matrix elements of Kuramoto et al. 7^ are related 
with our non-relativistic operators (|2.14p as: 



\if\m' 

l(/I^N>P 



47r^|(j;i||A^y||0+)p, 

4'^E E i(Jnii-M;^jiio+)p, 



J^ m=0,±l 



47r 



A = ^EO(J«ii-^ojiio 



+ \|2 



\{UMtj\\0+)\ 



+ \|21 



(A12) 



In Ref. [7J] are neglected the relativistic operators Mj, 
and M^j defined in (PJS)) . 



Appendix B: Sum Rule Approach 

We follow here the sum-rule approach developed by 
Kuramoto et al. [7J], and adapt it to our formalism. We 
start from Eqs. (|2.25p . and (|2.27p . and as in this work we 
assume that the wjtt dependence of the integrand can be 
ignored, fixing it at a representative value oTp" . The sum- 
mation over final nuclear states JJ^ then can be carried 
out by closure, and the ICS is 



ariE.) 



G^\Pll^F{Z + S,Ee) 



27r 



d{cosd)T 



SR 



(Bl) 



where the lepton energy is Eg ^ E^ — cuj^ , while the 
sum-rule matrix element reads: 



r^« = E {o+\oio^\o+)Ca 

a=0,O,±l 



25ft((O+|OjOo|O+)/:0o 



(B2) 



The operators Oa are given by p.6p . and the lepton 
traces by [H, Eq. (2.24)]. The matrix elements in (|B2|) 
arc proportional to A^(l — D), where iV^v is the number 
of neutrons (protons), contained in the target nucleus 
for the neutrino (anti-neutrino) reaction. The correla- 
tion functions D come from the Pauli-exclusion-effect, 
and depend on the type of the operator. One gets: 



V 



SR 



\ M 



2TfinC 



m'-iDo 



,(B3) 



with 



Tb = gUl - Ds) + {g^+g,,Y{l - Dl), 
To = t,{l ~ Ds) + {g^ -g,,f{l - Dl), 



Ti = (5A-ffw)'(l-^T), 
r_i = {g^ + g„f{l-DT), 



(B4) 



T(, 



00 



-3v5v(l - Ds) + (.gA +.9pi)(.9a - ffp2)(l - Dl). 



The correlation functions Ds,Dl and Ds were taken 
from t he S M calculation done by Bell, and Llewellyn 
Smith [ll6| with HO wave functions, and representing 
the nuclear ground state by a sin gle d eterminant wave 
function. The results for ^^C are Ill6l . Table ll): 



Ds 
Dt 



e-" [l + 0.14877^] , 

e-'' [0.704 + 0.14877 + 0.1487?^] , 



Dl = e"'' [0.704 + 0.29677 + 0.14877^] , (B5) 

where 77 = ^6^^^ ^0. 0558. 

As seen from (|2.26p . the factor |pf |i?f in (jBip behaves 
as {Eu — ujj^)^ , and therefore af^{Ei,) depends very crit- 
ically on the average value for the excitation energy aJp. 



Ei = E„-i^Ji,, |p£| = y'(£;^ -a;jj)2 -to2, 

2E^{Ei - Ip^I cos 6*) -ml+ oj]^, (B6) 



Appendix C: Extreme Relativistic Limit 

Using the present formalism the ERL, defined by the 
limit of the lepton velocity jp^j/i?^ — >■ 1, yields 

af^^(i?,) = y^^F{Z + S,E,) [ d(cos0)rj™^(«), 
^ 27r 7_i 

(CI) 
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with 



K = J2E^Ee{l -cos9)+ujj^, (C2) 



and 



2 cos^ - 



(J-||O0j(^)~^OojMI|O^ 



J2 i(jrj|o,„j('^)iio+)p 

m=±l 
k2 



where g^ = g^^ ~ 5pi- The muon capture transition rate 
reads 

H^m) = ^l'/'i5|'rA(c.j.), (D2) 

2 where (jjis is the muonic bound state wave function eval- 
uated at the origin, and 

+ 2|(j:||O_ij(ii;.)||0+)p], (D3) 



^ 1 ^^°^^2'^'^^''^^2 



with 



27n5 sin - 



e /F 



2 V K^ 



(C3) 



Appendix D: Muon Capture rate 

For the sake of completeness we also show the formula 
for the muon capture process within the present formal- 
ism. Here k ^ E^ = m^ — u;j» — AM — Eb, where E'^ 
is the binding energy of the muon in the IS* orbit, and 
instead of (|2.5p one has: 



K 



'2M' 



E^ 

^2M' 



9v — SvTjTj; <?A — 9a 

5w — [9v + gMJ-^TT', 5p — ffp^TT' 



(Dl) 



• For natural parity states , with tt = (—)'', i.e., J'^ 
0+,l-,2+,3-,---: 



, J- , i- , tv , 



and 



m^ — AE'coui — Eb . .V 
9v— ^- Mj , 

- (.9a + VjM^ij + g^.M'll, , (D4) 



• For unnatural parity states, with tt = i^)'^^^, i-G-, 
J'^ = 0-,l+,2-,3+,---: 

Obj-Oo.j = gAMl + {gA + gj,-gp)Moi, 

0-ij = -{gA + g„)Mt-"j-gvM]:{y (D5) 
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